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1 Introduction 

Let S be a smooth, closed Riemann surface of genus g, with n punctures and 
3g — 3 + n > 1. Let M_i denote the space of metrics of constant curvature 
— 1 on E. The group, Diff , of diffeomorphisms homotopic to the identity 
acts by pull back on M_ 1; and we can define Teichmiiller space T g n , to 

*This work was partially supported by a Nettie S. Autrey fellowship and NSF grants 
9971563 and 0139887. 
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be the quotient space M_i/ Dif /o(S). Note that the constant 3g — 3 + n 
is fundamental in Teichmiiller theory: it is the complex dimension of the 
Teichmiiller space T 9tn . 

The Weil-Petersson metric on Teichmiiller space has been heavily studied. 
It is defined by the L 2 inner product on the space of holomorphic quadratic 
differentials, the cotangent space of Teichmiiller space. This metric has many 
curious properties; for example, it is a Riemannian metric with negative 
sectional curvature EH]), but it is incomplete (jjj [2E])- Also this metric 
is Kahler (|2J), and there is an upper bound 2 7i-(g-i) ^ or ^ ne holomorphic 
sectional curvature and Ricci curvature (|2H|)> a result conjectured by Royden 
(|19j). However, there is no known negative upper bound for the sectional 
curvatures. 

In this paper, we investigate the asymptotic behavior of the sectional cur- 
vatures of the Weil-Petersson metric on Teichmiiller space. Our method is to 
investigate harmonic maps from a nearly noded surface to nearby hyperbolic 
structures. One of the difficulties in estimating Weil-Petersson curvatures 
is working with the operator D = — 2 (A — 2) _1 , which appears in Tromba- 
Wolpert's curvature formula. Our approach is to study the Hopf differentials 
associated to harmonic maps and the analytic formulas resulting from the 
harmonicity of the maps. There is a natural connection between the operator 
D and the local variations of the energy of a harmonic map between surfaces 
(|25j). From this connection, we estimate the solutions to some ordinary 
differential equations to derive our curvature estimates. Our estimates imply 
that even though the sectional curvatures are negative, they are not staying 
away from zero. More specifically, we prove our main theorem: 

Theorem 1. If the complex dimension of Teichmiiller space T is greater than 
1, then the sectional curvature of the Weil-Petersson metric is not pinched 
from above by any negative constant, i.e, there is no negative upper bound 
for the Weil-Petersson sectional curvature. 

We remark that this asymptotic flatness results from pinching two inde- 
pendent closed geodesies on the surface. We will discuss more asymptotic 
flatness in an upcoming paper. The proof of the main theorem provides a 
quantitative bound for the sectional curvature, namely, 

Theorem 2. If the complex dimension of Teichmiiller space T is greater 
than 1, and let I be the length of shortest geodesic along a path to a boundary 
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point in Teichmiiller space, then there exists a family of tangent planes with 
Weil-Petersson sectional curvature of the order 0(1) . 

Recently Brock and Farb ([5 ) showed that the Weil-Petersson metric on 
T is Gromov hyperbolic if and only if dim c {T) < 2. They pointed out that 
if the Weil-Petersson metric on T had curvature pinched from above by a 
negative constant, the Gromov hyperbolicity in the case of the surface being 
doubly-punctured torus or 5-punctured sphere would be an immediate conse- 
quence of the comparison theorems. In the end of their paper, they asked the 
following question: if mi(£) is homeomorphic to a doubly-punctured torus or 
5-punctured sphere, are the sectional curvatures of the Weil-Petersson metric 
bounded away from zero ( 5J)? Our paper gives a negative answer to this 
question. 

Naturally associated to a harmonic map w : (T,,a\dz\ 2 ) — > {Y,, p\dw\ 2 ) is 
a quadratic differential $(cr, p)dz 2 , which is holomorphic with respect to the 
conformal structure of a. This association of a quadratic differential to a 
conformal structure then defines a map $ : T — ► QD(E) from Teichmiiller 
space T to the space of holomorphic quadratic differentials QD{Ti). This 
map is known to be a homeomorphism (|25j). 

As an important computational tool in geometry of Teichmiiller theory, 
the method of harmonic maps has been studied by many people. In particu- 
lar, the second variation of the energy of the harmonic map w = w(o~, p) with 
respect to the domain structure a (or image structure p) at a = p yields the 
Weil-Petersson metric on T ([2H], [23), and one can also establish Tromba- 
Wolpert's curvature formula of the Weil-Petersson metric from this method 

([131, see also 

The moduli space of Riemann surfaces admits a compactification, known 
as the Deligne-Mumford compactification ([E]), and any element of the com- 
pactification divisor can be thought of as a Riemann surface with nodes, a 
connected complex space where points have neighborhoods complex isomor- 
phic to either {\z\ < e} (regular points) or {zw = 0; \z\, \w\ < e} (nodes). 
We can think of noded surfaces arising as elements of the compactification 
divisor through a pinching process: fix a family of simple closed curves on E 
such that each component of the complement of the curves has negative Eu- 
ler characteristic. The noded surface is topologically the result of identifying 
each curve to a point, the node (jl]). 

While writing this paper, the author learned that Scott Wolpert f[3Ti]) 
has obtained results related to ours. 
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The organization of this paper is as follows. In section 2 we give the 
necessary background, define our terms and introduce the notations. Section 
3 is devoted to proving our main theorem. The discussion of this purpose 
is broken into subsections: in §3.1, we collect the local variational formulas 
associated with the harmonic maps, hence connect the local variations of the 
energy of the harmonic map to the operator D = — 2(A — 2)~ x ; in §3.2 we 
describe the so called "model case" of the problem; namely, we pinch two 
core geodesies of two cylinders into two points, and study the asymptotic 
behavior of the harmonic maps between cylinders. In the model case of this 
problem, the surface is a pair of cylinders; in §3.3, we describe a family of 
asymptotically flat 2-dimensional subspaces. We also establish two ordinary 
differential equations, and relate the solutions of the equations to the operator 
D = — 2(A — 2) _1 , which leads us to estimate the curvature. We will establish 
the estimates of terms in the curvature formula when the surface is a pair 
of cylinders; in §3.4 we construct families of maps which have small tension, 
and are close to the harmonic maps resulting from pinching process in §3.3; 
finally in §3.5 we prove our main theorem using the estimates in §3.3 and the 
construction in §3.4. In section 4, we comment on the case when the surface 
has finitely many punctures. 

Acknowlegement: This paper contains part of the author's doctoral 
dissertation at Rice University. It is his great pleasure to express his deepest 
gratitude to his thesis advisor, Michael Wolf, for suggesting this problem, 
many stimulating conversations and for his patience. He also wants to thank 
Professor Frank Jones for his generous help. 

2 Notations and Background 

Recall that £ is a fixed, oriented, smooth surface of genus g > 1, and n > 
punctures where 3g — 3 + n > 1. We denote hyperbolic metrics on £ by 
a\dz\ 2 and p\dw\ 2 , where z and w are conformal coordinates on £. By the 
uniformization theorem, the set of all similarly oriented hyperbolic structures 
M_i can be identified with the set of all conformal structures on £ with the 
given orientation. Equivalently, this is the same as the set of all complex 
structures on £ with the given orientation. And Teichmuller space T is 
defined to be 

T= M.jDiffoiE) 
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It is well-known that Teichmiiller space T has a complex structure (j2]) 
and the cotangent space at a point E G T is the space of holomorphic 
quadratic differentials Qdz 2 on E. On E, there is a natural pairing of 
quadratic differentials and Beltrami differentials fJ>(z)jr- given by 

</i, $ >= Re f^p(z)&(z)dzdz 

The tangent space at E is the space of Beltrami differentials modulo the 
trivial ones which are ones such that < p, $ >= 0. The Weil-Petersson 
metric on T is obtained by duality from the L 2 -inner product on QD(E) 

< (j), if) >= L ^-dzdz 

where adzdz is the hyperbolic metric on E. 

For a Lipschitz map w : (E, <r|<i2:| 2 ) — > (E, p|rfw| 2 ), we define the energy 
density of w at a point to be 

e(w) = ^KI 2 + ^K"| 2 

and the total energy E(w; a, p) = L e(w; cr, p)o~dzdz. 

A critical point of E(w;o~,p) is called a harmonic map; it satisfies the 
Euler-Lagrange equation, namely, 

Wzz + e fvo z w 2 = 0. 

The Euler-Lagrange equation for the energy is the condition for the van- 
ishing of the tension, which is, in local coordinates, 

t(w) = Aw"' + N Tlpwfwf = 

It is well-known [H] [H] |22] 1201) that given a, p, there exists a unique 
harmonic map w : (E, cr) — > (S,p) homotopic to the identity of E, and this 
map is in fact a diffeomorphism. 

Thus we have a holomorphic quadratic differential (also called a Hopf 
differential) §dz 2 = pw z w z dz 2 , and evidently 

$ = w is conformal ^ a = p. 

where a = p means that (E, a) and (E, p) are the same point in Teichmiiller 
space T . This describes $ as a well-defined map $ : T — » QD(11) from 
Teichmiiller space T to the space of Hopf differentials QD{Yi). This map 
is a homeomorphism (|25|). Also note that the map w can be extended to 
surfaces with finitely many punctures. 

We define two auxiliary functions as following: 
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n 



= H{<r,p) = 



Hence the energy density is e = Tt + C Many aspects connected with a 
harmonic map between Riemann surfaces can be written in terms of 7i and 
£ (and $). We denote on (£,cx|dz| 2 ) 

A = tP^Ktp) = -^^logp,K(a) = -^^loga, 

a azaz' \ r ' p awow • 3r ' v / a azaz & ' 

where A is the Laplacian, and K(p), K(a) are curvatures of the metrics p 
and a, respectively. 

The Euler-Lagrange equation gives 

AlogH = -2K(p)H + 2K(p)C + 2K(a). 

When we restrict ourselves to the situation K(a) = K{p) = — 1, we will 
have the following facts ([25]): 

• H > 0; 

• The Beltrami differential a = — = 

™ w z art ' 

• AlogH = 2H-2£-2. 

3 Weil-Petersson Sectional Curvature 

In this section, we will prove the main theorem of this paper, namely, 

Theorem 1. If the complex dimension of Teichmiiller space T is greater than 
1, then the sectional curvature of the Weil-Petersson metric is not pinched 
from above by any negative constant, i.e, there is no negative upper bound 
for the Weil-Petersson sectional curvature. 

Combined with the fact that there is no lower bound for the sectional 
curvature ([21]), We have 

Corollary 3. There are no negative bounds for the sectional curvature of the 
Weil-Petersson metric when dimc(T) > 1. 



6 



The rest of this paper will be devoted to proving our main theorem. 
The discussion will be broken into 5 subsections. In §3.1, we collect some 
variational formulas we will need in the rest of discussion; in §3.2, we mostly 
describe the "model case", i.e. we consider a family of pertubations of the 
identity map between a pair of cylinders; in §3.3, we describe a family of 
subplanes of the tangent space of Teichmiiller space T gi2 , and estimate terms 
in the curvature formula when the surface is a pair of cylinders; in §3.4, we 
construct C 2,a maps between surfaces and show that the constructed maps 
have small tension and are close to the harmonic maps we get from the 
pinching process; in §3.5, we prove the main theorem. 

3.1 Local Variation 

We consider a family of harmonic maps w(t) for t small, where w(0) = id, 
the identity map. Denote by $(£) the family of Hopf differentials determined 
by w(t). We can rewrite the equation AlogH = 2H — 2C — 2 as 

AlogH(t) = 2H(t) - ^ggf - 2 

For this equation, we see that the maximum principle will force all the odd 
order t-derivatives of TC(t) to vanish, since the above equation only depends 
on the modulus of $(t) and not on its argument (j25*j), and T~C{t) is real- 
analytic in t (j2E])- 

Wolf computed the i-derivative of various geometric quantities associated 
with the harmonic maps, and we collect these local variational formulas into 
our next lemma: 

Lemma 4. f}25f ) For the above notations, we have 

• Kit) > I, and Hit) = 1 4» t = 0; 

• U{t) =d/dt a \ H(t) =0; 

• ii = d/dt a \ ^{t) = ®Ja; 

• ii{t) = Q^\aH{t) = D^r 1 , where D = -2(A - 2)" 1 . Evidently D 
is a self-adjoint compact integral operator with a positive kernel, and it 
is the identity on constant functions. 

With this lemma, we obtain a partial differential equation about TC(t): 

(A-2)(W(f)) = -25£* 
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3.2 Harmonic Maps Between Cylinders 

For the sake of simplicity of exposition, we assume that our surface has no 
punctures. We will comment on punctured case in section 4. 

Before we jump into the proof of the main theorem, in this subsection, 
we consider the asymptotics of harmonic maps when our surface is a pair of 
cylinders. We consider the asymptotics of a family of pertubations of the 
identity map between two cylinders, and we call this is the "model case" . In 
the model case, the surface is a pair of cylinders M and Mi. In particular, 
consider the boundary value problem of harmonically mapping the cylinder 

M = [l^sin-^l), nl- 1 - l-\sin- l {l)} x [0, 1] 

with boundary identification 

where the hyperbolic length element on M is lcsc(lx) \dz\, to the cylinder 

N = [L^sin-^L), tiL- 1 - L^sirT^L)] x [0, 1] 
with boundary identification 

where the hyperbolic length element on iV is Lcsc(Lu) \dw\. Here I and L are 
the lengths of the simple closed core geodesies in the corresponding cylinders. 
Say w = u + iv is this harmonic map between cylinders M and N, where 

u(l, L; x, y) = u(l, L; x),v(x, y) = y. 

The Euler-Lagrange equation becomes 

u" = Lcot(Lu)(u' 2 - 1). 

with boundary conditions u( sm l — ) = sm L — and m(^) = Note that 
both M and iV admit an anti-isometric reflection about the curves {^} x [0, 1] 
and {^} x [0,1]. 

Since the quadratic differential $ = pw z w z = \L 2 csc 2 (Lu){u' 2 — 1) is 
holomorphic in Mo, we have 

= §(p 2 w z w z ) = f x {\L 2 csc 2 {Lu){u' 2 - 1)) 
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Therefore L 2 csc 2 (Lu)(u' 2 — 1) = c (l,L), where c (l,L) is independent of 
x, and Co(/, Z) = since u(x) is the identity map when L = I. 
So we have 

u' = a/1 + co(/, L)L~ 2 sin 2 (Lu) 

with boundary conditions u( sm -, — ) = — L — and ^(f^) = So the 

solution to the Euler-Lagrange equation can be derived from the equation 

f" lr . d " = z-i- 1 sm- 1 m 

JL ism 1 L ^i +CQ (i,L)L 2 sin 2 {Lv) K ' 

with co(LL) chosen such that f r 2 fi . _ ir , = = = h — l~ 1 sin~ 1 l. 

uv ' ' JL sm L y/l+co(l,L)L 2 sin 2 (Lv) 21 

It is not hard to show that when I — > 0, the solution u(l, L; x) converges 
to a solution u(L;x) to the "noded" problem, i.e., M = [l,+oo), where we 
require u(L] 1) = and \im x ^ +OCl u(L; x) = j^. 

This "noded" problem has the explicit solution as following (|26.) 

i r 2L(l-x) 

u(L;x) = L^sin- \^ e2L(1 _ x) } 

- L_r (l+L) e 



with the holomorphic energy 

H (L;x) 



„ „ 1 1 ./ BES e&q-*) 

L 2 z 2 r 1+ V (i+i-) e i 

4 [ i_,/iiEir e i(i-x) J 

V (i+i) e 



3.3 Asymptotically Flat Subplanes 

Consider the surface S which is developing two nodes, i.e., we are pinching 
two nonhomotopic closed geodesies 7 and 71 on S to two points, say p and 
Pi. We denote M and Mi their pinching neighborhoods, i.e., two cylinders 
described as M above centered at 7 and 71, respectively. 

We define M(Iq, l\) be the surface with two of the Fenchel-Nielsen coordi- 
nates, namely, the hyperbolic lengths of 70 and 71, are Iq and l±, respectively. 
When we set the length of these two geodesies equal to I simultaneously, we 
will have a point M(l) = M(l, I) in Teichmuller space T g . Note that as I 
tends to zero, the surface is developing two nodes. At this point M(l), we 
look at two directions. First, we fix 71 in Mi having length I, and pinch 70 
in M into length L = Lit), where L(0) = I. So the t-derivative of /x (t), the 
Beltrami differential of the resulting harmonic map, at t — represents a 
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tangent vector, fi , of Teichmiiller space T g at M(l); we denote the resulting 
harmonic map by Wo(t) : M(l, I) — > M(L(t), I). Then we fix 70 in M having 
length I, and pinch 71 in Mi into length L = L(t), so the t-derivative of 
/ii(t) at t — represents another tangent vector, fii, at M(Z); we denote the 
resulting harmonic map by W\{t) : M(l, I) — > M(l; L(t)). These two tangent 
vectors /i and fii will span a 2-dimensional subspace of the tangent space 
Tm(i)% to T g , hence we obtain a family, Qi, of 2-dimensional subspaces of 
the tangent space of T g . 

Proposition 5. The Weil-Petersson sectional curvatures of Qi tend to zero 
as I tends to zero. 

It is immediate that this proposition implies our main theorem. 
To begin the proof of the proposition, we look at the curvature tensor of 
the Weil-Petersson metric, which is given by f|29p 

where dA is the area element. Then the curvature of Q is R/H, where ([29 ) 
R — Rqioi — Rqiw ~~ -R1001 + -R1010 

and 

II = 4 < fi , fi >< > — 2| < /i , fix > | 2 — 2Re(< {iq, fii >) 2 

= 4 < /i ,Ao >< > -4 I < A*o,Ai > I 2 

In the rest of this subsection, we will estimate terms in the above curva- 
ture formula in the model case, i.e., the surface is a pair of cylinders M and 
Mi. 

We denote <f)o(t) as the Hopf differential corresponding to the cylinder 
map Wo(t) in Mq, and <pi(t) as the Hopf differential corresponding to the 
cylinder map ui\{t) in Mi. Here Wo(t) : M (7) — > M (L(t)) and u>i(t) : 
M(0 ~ * M (£(£)) are harmonic maps described as w = u(x) + iy in the 
last subsection. In Mi, we still use O to denote the Hopf differential corre- 
sponding to harmonic map Wo(t), while in Mq, we still use <pi as the Hopf 
differential corresponding to harmonic map W\(t). We also denote /xo and 
fii as the corresponding Beltrami differentials to 0o and 0i- 

We denote a = a(l) = / _1 sm _1 (Z), and b = b(l) = 7r/ _1 — / _1 sm _1 (Z). 
And in this paper, A ~ i? means v4/C<-B<Cv4 for some constant C > 0. 
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So in M , as in §3.2, we choose L(t) so that ^\ t=0 co(t) = -^\l=iCo(1,L) = 
4. Here we recall that co(l,L(t)) = L 2 csc 2 (Lu)(u' 2 — 1). We notice that do 
is never zero for all / > 0, otherwise we would have w z = as = aw z . 
Hence w is a constant map by rotational invariance of the map. Thus we 
have 0o = f t \t=o(\co(t)) = 1 in M . 

Also, we see that |0i||m o = f° r x £ [ a , b], where ((x, I) satisfies that 
C(x,l) < C\x~ A for x E [a,n/2l], and COM) < Ci(n/l -a:)" 4 for x E [vr/2/,6], 
and ((x, 0) decays exponentially in [1, +oo]. Here C\ = C\(l) is positive and 
bounded as / — > 0. To see this, notice that 0i is holomorphic and \(f>i\ is 
positive, so log\(j>i\ is harmonic in the cylinder M . Hence we can express 
log\4>\\ in a Fourier series T 1 a n {x)exp{— iny), and we compute = A/og|0i| = 
S(a^ — n 2 a n )exp(—iny) . We will see, in §3.4, that 0i is close to in M Q . 
Hence we conclude the properties ( has. Similarly, we have 0i|mi = 1 and 
I0o||mi = C(x,l) for x E [a,b]. 

Note that these estimates imply, informally, that most of the mass of \4>o\ 
resides in the thin part associated to 70, and most of the mass of \<f>i\ resides 
in the thin part associated to 71. 

Now in Mo, the corresponding Beltrami differential is 

Ao = = io/v 

and |/i | 2 |M„ = |0oH 2 |m o = / _4 sm 4 (/a;). 
Also, in that same neighborhood M , 

/tl = 0i/(7 

Hence 

IAi| 2 |m = |0i/o"| 2 |m o = l~ i sin i (lx)C 2 (x,l). 
Similarly, in M ± , we have 

IAo| 2 |mi = l~ 4 sin' l (lx)( 2 (x,l) 
|/*i| 2 |mi = l~ A sin A {lx) 

Lemma 6. 1/11 = 0(Z 3 ). 

We recall that n = 4 < /io,Ao >< > — 4| < /ioiAi > | 2 5 and 

compute the asymptotics in / of each term. Using |/io| 2 |M — l~ A sin A (lx), 
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and noting that a = a(l) = I l sin and b = b(l) — irl 1 — I 1 sin 

we have that 



<Ao,Ao>|m = / \fi \ 2 crdxdy 
J M 

= / / \fio\ 2 °~dxdy 

JO J a 



1 r b 

-2-2 




/ sin Ixdxdy 



o J a 



= ![/-3_ r 3 sm -l (/) 

2 v 1 

7-3 

And using |/ii| 2 |M = l~ A sin 4 (lx)( 2 (x, I), and Z) < Cia;~ 4 forx G [a, 7r/2Z], 
we have 



</ii,/ii>| Mo = / |/ii| 2 <7<Zzdy 



< 



M 

2 / / r 4 sin A (lx)( 2 (x,l)adxdy 

JO J a 

pi pit/ 21 

2Cl / / r 2 sin 2 (lx)x- 8 dxdy 

Jo Ja 



= 0(1) 

Also < fiQifii >= Jj.fioP'idA, hence, 

< fio, > \m = / fiofiivdxdy 

J Ma 



< Ci I 2 sin 2 (lx)x A dxdy 

J Mo 



Similarly, 



'Mo 

0(1) 



< /i , Ao > |mi = 0(1) 

< > | Mi ~ l~ 3 

< /i ,/il > | Mi = 0(1) 
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Note that LT > (4 < /io,Ao >< > — 4(< /ioiAi >)) 2 

which completes the proof of Lemma 6. 
From Lemma 6, we have 



Mo 



\R\/n = o(\R\/(r 3 )) = o(\R\i 3 ) (i) 

Now we are left to estimate \R\. We have 
Lemma 7. \R\ <4 J E D(\fi \ 2 )\p,i\ 2 adxdy 

Note that D = — 2(A — 2) _1 is self-adjoint, hence we have 
J s Z)(|/io| 2 )|Ai| 2 cr<ir<% = J^D(\fii\ 2 )\fi \ 2 adxdy 

Therefore, 

R — -Roioi ~~ -Roiio — -Riooi + -Rioio 

= 2 D(fj, j2i)fj, fiiadxdy + 2 / D(tiijj,o)iiijio<jdxdy 



D(\/l \ 2 )\jli\ 2 adxdy — ^ D(|/ii| 2 )|/io| 2 <TG?x<i?/ 
D(ji p ll )fi 1 p Q adxdy- / Difi^/lo^adxdy 



= 2 D(fioiii)iiifi adxdy — 2 / Z}(|/i | 2 )|/ii| 2 a<i£cfo/ 
7s is 

The last equality follows from that here /io and /ii are real functions. Now 

from lemma 4.3 of [HJ, we have |£>(/i /ii)| < |#(|/i | 2 )| 1/2 |£(|/ii| 2 )| 1/2 . Then 

an application of Holder inequality shows that 

| / DitiofiJfi^odAl < [ iDd/ioHr^Pd/iiHr^/ii/io^ 

D(|/io| 2 )|/ii| 2 cL4)s(/ Dd^l 2 )^! 2 ^)* 



D(|/io| 2 )|/ii| 2 cr<ix(i?/ 



So we will have 



\R\ <4 / D(|/i | 2 )|/ii| 2 flrdxdy (2) 

which completes the proof of Lemma 7. 

Therefore together with (1) and (2), to show the proposition, it is enough 
to show that J E D(\ii \ 2 )\iii\ 2 adxdy = o(l~ 3 ). In fact, 
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Lemma 8. J T; D(\ii \ 2 )\fi 1 \ 2 adxdy = O(l 2 ) 

Firstly, from Lemma 4, we have 

D(|/i | 2 ) = -2(A-2)- 1 ^ 

We recall in M , the Hopf differential O is corresponding to the cylinder 
map w (t) : (M ,a) — > (M ,p(t)), and \<j> \ = 1. As in section 2 and §3.1, 
we write 7i = | | 2 , therefore we can write D(\jl \ 2 ) = H. Then in the 

cylinder M , we have 

(A - 2)H = -2^1! = -2r A sm\lx) (3) 
a 2 

A maximum principle argument implies that 7i is positive. This 7i con- 
verges to the holomorphic energy H in the "noded" problem (of §3.2) when 
x is fixed but sufficiently large in [a,ir/2l]. This convergence guarantees 
that 7i is bounded on the compacta in [a,b]. Thus we can assume that 
Ai{l) = H(a) = liil^sin'H) = 0(1) > 0. Then H(x) solves the following 
differential equation 

(l- 2 sin 2 {lx))H" -2H= -2l- 4 sin\lx) (4) 

with the conditions 

Hil^sin-H) = A 1 (l),H'(n/2l) = 

Recall from §3.1 that all the odd order t-derivatives of H(t) vanish. Also 
notice that 

J{x) = 

is a particular solution to equation (4). Hence we can check, by the method 
of reduction of the solutions, the general solution to equation (4) with the 
assigned conditons has the form 

H(l; x) = J(x) + A 2 cot(lx) + A 3 (l - lxcot(lx)) 

where coefficients A 2 = A 2 (l) and A 3 = A 3 (l) are constants independent of x 
and we can check, by substituting the solution into the assigned conditions, 
that they satisfy that 
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A 2 = f A 3 = 0(1-') 
Noticing that ((x,l) < C\X~ A in [a,ir/2l], we compute the following 



7i{x)\^ii\ adxdy = 7i{x){l~ 2 sin 2 {lx))C > 2 {x,l)dxdy 

JO J a 

rl ptt/21 

< 2C 2 ( / / J(x)r 2 sin 2 (lx)x- 8 dxdy 

J0 J a 



+ A 2 cot(lx)l~ sin (lx)x dxdy 

Jo J a 

+ / / A 3 (l - lxcot(lx))l~ 2 sin 2 (lx)x~ 8 dxdy) 

Jo J a 

-2\ , na-^ 



= o(r 2 ) + o(r 2 ) + o(r 2 ) 
= o(i- 2 ) 



Therefore 



2 f b - 2 

-D(|/i | 2 )|/ii| adxdy = / H\fi\\ adx 

j J a 

= 0(r 2 ) (5) 

Now let us look at the term J M D(\fio\ 2 )\fii\ 2 adxdy. Note that in Mi, 
which we identify with [a, b] x [0, 1], we have D(\fi \ 2 ) = H, here fio(t) and H 
come from the harmonic map W (t) : M(l,l) — * M(L(t),l), where | /io 1 1 iWi = 
l~ 2 sin 2 (lx)((x,l), and H(t) solves 

(r 2 sin 2 {lx))H" -2H = -2r 4 sin\lx)( 2 {x, I) (6) 

with the conditions 

Hil-Ksin-H) = Bi(l), H(tt/1 - l^sin-H) = B 2 (l). 

Here B ± (l) and B 2 (l) are positive and bounded as / — > 0, since Ti converges to 
the holomorphic energy in the "noded" problem (of §3.2, when M x = [1, oo)). 
We recall that |0i||m o = ((x,l), where ((x,l) < Cix~ A for x e [a, 7r/2Z] and 
t(x, I) < C^n /I - x)~ 4 for x G [n/2l, 6]. 
Consider the equation 

(r 2 sin 2 (lx))Y" -2Y = (7) 
with the boundary conditions that satisfy 
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Y^sirrH) = 0(1) and Y(n/l - l^siwH) = 0(1), as Z -> 0. 

We claim that there exists some h = 0(1) such that H — h is a supersolution 
to (7) for x G [/~ 1/4 , 6-Z -1 / 4 ]. To see this, we notice that 2\<j) \ 2 \ Ml l~ 4 sin 4 (lx) 
decays rapidly as x — > 7r/2Z for small Z. So there is a positive number ft, = 
0(1) such that 2|0 o | 2 | Ml Z~ 4 sm 4 (Z:r) < 2h for x G [l- l/ \b - Now for 

x G [Z -1 / 4 , 6] we have 

(r 2 sin 2 (lx))(H - h)" - 2(H -h) = 2h- 2\^\ 2 \ Ml l~ i sin\lx) > (8) 

Notice that for any constant A, we have that if Y(x) solves the equation 
(7), then so does \Y(x). So up to multiplying by a bounded constant, we 
have Y\du 1 > (H — h)\gM l - Hence Ti — h is a subsolution to (7) for x G 
[Z -1 / 4 , b — Z -1 / 4 ]. We can check the solutions to (7) have the form of 

Y(l; x) = B 3 cot(lx) + B 4 (l - lxcot(lx)) 

where constants B 3 = B 3 (l) and _B 4 = B 4 (l) satisfy, from the boundary 
conditions for the equation (7), that 

B 3 = 0(1) and B 4 = 0(1) 

Therefore in [Z" 1 / 4 , b - Z" 1 / 4 ], We have H < h + Y(x). Now, 

/ Y(x)\fli\ 2 adxdy — Y(x)(l~ 2 sin 2 (lx))dxdy 

JMx Jo J a 

B 3 cot(lx) (l~ 2 sin 2 (lx))dxdy 



i 

J a 



+ B 4 (l — lxcot(lx))(l 2 sin 2 (lx))dxdy 

JO J a 



= o(r 2 ) + o(i- 2 ) 
= o(r 2 ) 

Also for x G [a, Z -1 / 4 ] U [b — Z -1 / 4 ,6], we apply maximum principle to 
equation (6) and find that 

H < max(l- 4 sin 4 (lx )( 2 (x , 1), sup(H\ matl -i/4 ]u[b _ l -i/4 M)xm) )) 

for some x G [a, Z -1 / 4 ]. Hence using the properties of ((x, I) in [a, b], a direct 
computation shows that 
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SiS^n(x)\^ 1 \ 2 dA = o(i-^) = {i- 1 ) 

With this and using H < h + Y(x) for x G [Z -1 ^ 4 , b — / _1 ^ 4 ], we have 

-1 r b 



I D(\ii Q \ 2 )\^\ 2 dA = f f 1i(x)\^\ 2 dA 

J Mi JO Ja 

_ l 

pi pb—l 3 pi pi t 

/ / H(x)\fi 1 \ 2 dA+ / / 7i(aO|fii| 2 cL4 

JO J/- 1 /4 Jo Ja 



+ [ [ H(x)\^\ 2 dA 
Jo Jb-1- 1 / 4 



< [ [ (Y(x) + h)\fi 1 \ 2 adxdy + o(r 1 ) 

Jo Ja 

= o(r 2 ) + o(r 2 ) + o(r 1 ) = o(r 2 ). (9) 

Now combine the estimates of Lemma 7, (9), (10), we will have 
\R\ < A J D(\^\ 2 )\^\ 2 adxdy 

= A(f D(\ f i \ 2 )\ f i l \ 2 dA+ [ D(\ii \ 2 )\»i\ 2 dA 

J M J Mi 

= o(i- 2 ) (10) 

Now with Lemma 6 we have 

\R\/Il = 0(l~ 2 )l 3 = 0(1) 
namely, |-R|/IT tends to zero as we pinch the geodesies. 

3.4 Constructed Maps 

In last subsection, we estimated the terms in the curvature formula in the 
model case, i.e., the surface is a pair of cylinders. Essentially, in each pinching 
neighborhood Mj, we used the cylinder map Wi instead of the harmonic map 
Wi during the computation, for i — 0,1. Now we are in the general setting, 
i.e., the surface S develops two nodes. In this subsection, we will construct 
a family of maps Gj to approximate the harmonic map Wi, and the essential 



17 



parts of this family are the identity map of the surface restricted in the non- 
cylindral part and the cylinder map in each pinching neighborhood. We will 
also show that this constructed family Gi is reasonably close to the harmonic 
maps Wi, fori = 0,1; hence we can use the estimates we obtained in the 
previous subsection to the general situation. 

We recall some of the notations from previous subsections. We still set M 
and Mi to be the pinching neighborhoods of the nodes po and pi, respectively. 
Also Wo(t) is the harmonic map corresponding to fixing 71 in Mi having 
length /, and pinching 70 in Mo into length L = L(t), where L(0) = I. 
Similarly W\{t) is the harmonic map corresponding to fixing 70 in M into 
length /, and pinching 7! in Mi having length L = L(t). Let w (t) and Wi(t) 
be cylinder maps in the model case, we want to show that W (t) is close to 
wo(t) in Mo and is close to identity map outside of M> 

We denote subsets E = {p E E : dist(p,dM ) > 1}, and Si = {p E 
E : dist(p,dMi) > 1}. Define the 1-tube of dM as B(dM ,l) = {p E 
E : dist(p,dM ) < 1}, and the 1-tube of dM x as B{dM u \) = {p E E : 
dist{p, dMi) < 1}. We can construct a C 2 ' a map G : E — > E such that 

f Wo(*)(p), pGM nE 
Go(p) = < P, pe (So\M ) 
[ y t (p), P eB(9M ,l) 

where gt(p) in B(dM , 1) is constructed so that it satisfies: 

• 9tip) = P for p E 9(E \(M U S(9M ,1))), and g t {p) = w (t)(p) for 

pea(M nE ) ; 

• ^ is the identity map when t — 0; 

• ^ is smooth and the tension of <7 t is of the order 0(t). 

We note that G consists of three parts. It is the cylinder map of M in 
Mo H So, the identity map in Eo\Mo, and a smooth map in the intersection 
region B(dMo, 1). Among three parts of the constructed map Go(t), two of 
them, the identity map and the cylinder map, are harmonic, hence have zero 
tension; thus the tension of G (t) is concentrated in B(dM , 1). From §3.2, 
for the cylinder map w = u(x) + iy, we have u' = a/1 + c (t)L~ 2 sin 2 Lu, 
where c (0) = 0, c (0) = 4,. Hence for x E [I' 1 sin' 1 (I) , I" 1 sin" 1 (I) + 1], 

w , z (x, y) = \{u\x) + 1) = |((1 + 0(l)cb(*))5 + 1) = 1 + 0(l)t + Oit 2 ) 
\w , z (x,y)-l\ = O(t) -> 0, (t -> 0) 
\w 0tZ z(x,y)\ = \\u"(x)\ = 0(\Lcot(Lu)(u' 2 - 1)|) = Q(t) 
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Thus we can require that \gt, z \ < C 2 t and \gt,zz\ < ^t, and the constant 
G% = C 2 (t, I) is bounded in both t and / since the coefficient of t for c (t) is 
bounded for small t and small I. With the local formula of the tension in 
section 2, we have r(Go(t)), the tension of Go(t) is of the order 0{t). 

Similarly, we have a C 2 ' a map G\ : S —>■ S. Note that these constructed 
maps Go and G\ are not necessarily harmonic. 

Now we are about to compare the constructed family Go(£) and the family 
of the harmonic maps Wo(t). To do this, we consider the function Qq = 
cosh(dist(W ,G )) - 1. 

Lemma 9. dist(Wo, Go) < C^t in B(dM ,l), where the constant G3 = 
Cz{t, I) is bounded for small t and I. 

First, we want to show that Qq is a C 2 function. Notice that both the 
harmonic map W (t) and the constructed map G (t) are the identity map 
when t = 0, and both families vary smoothly without changing homotopy 
type in t for sufficiently small \t\ (|E])- For all I > 0, and for any e > 
0, there exists a 5 such that for \t\ < 5, we have |Wo(t) — Wo(0)| < | 
and \Go(t) — Go (0)| < |. Therefore the triangular inequality implies that 
|Wo(t) — Go(t)| < e. Since I is positive, the Collar Theorem (0) implies that 
the surface has positive injectivity radius r bounded below, and we choose 
our e « r, so Q is well defined and smooth. 

We follow an argument in [T2]. For any unit v G T 1 (5(9Af , 1)), the map 
G satifies the inequality |||dG (u)|| — 1| = 0(t), hence \dG (v)\ 2 > 1 — e 
where Eq = 0(t), then we find that for any x e S, 

AQ > min{\dG (v)\ 2 : dG (v)±. lx }Q 

- <T{G ),X7d(;W )\ Go{x) > sinh(dist(W ,G )) (11) 

where 7^ is the geodesic joining Gq(x) to Wq{x) with initial tangent vector 
— v d(», W )\g ( x ) and terminal tangent vector \/d(G (x), •)|wo(«)- ^ ^o(^) 
does not coincide with Wo(t) on B(dM , 1), we must have all maxima of Qo(t) 
on the interior of B(dM , 1), at any such maximum, we apply the inequality 
\dG (v)\ 2 > 1 - e to (12) to find 

0>AQ > (l-e )Q -T(G )(sinh(dist(W Q ,G )) 
so that at a maximum of Qq, we have 

^ r(Go)sinhdist(Wo,Go) 

^° - (T^) 
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We notice that Q is of the order dist 2 (W Q , G ) and sinhdist(Wo, G ) is 
of the order dist(W , G ), this implies that dist(W , G ) is of the order O(t) 
in B(dMo, 1), which completes the proof of Lemma 8. 

Remark 10. Lemma 8 implies that Qo(t) is of the order 0(t 2 ) in B(dM , 1). 

Note that B(dM , 1) contains the boundary of the cylinder M nE , which 
we identify with [a + 1,6 — 1] x [0, 1], where, again, a = a(l) = / _1 sm _1 (/), 
and 6 = b(l) = nl^ 1 — l~ l sin~ l {l). While in the cylinder Mo H So, we have 
the inequality 

AQ > (l-e )Q -T(G )(sinh(dist(W ,G )) 

= (1 - e Q )Q - r(G )(tanh(dist(W , G ))(l + Q ) 
= (1 - so - r(G ))Q - r(G )(tanh(dtst(W , G )) 
> 1/2Q -C 4 t 2 

where the constant C4 is bounded for small t and /. Therefore we find that 
Qo(z,t) decays rapidly in z = (x,y) for x close enough to tt/21. Hence we 
can assume that dist(W , G ) is at most of order C't in [a + l,b — 1], here 
C = C'(x,l) is no greater that C 5 x~ 2 for x G [a + 1,tc/21], and no greater 
than C^tt/I — x)~ 2 for x G [ir/2l,b — 1], where C5 is bounded for small t 
and I. Both maps Wo and Go are harmonic in Mo D So, so they are also 
C 1 close (jHl), i.e. we have \W 0j2 — Go,z\ < C 5 x~ 2 t for small t and I, when 
x G [a + l,7r/2/]. Thus we see that |Wq,z — Go^l = C 5 x~ 2 , for x G [a + l,7r/2Z]. 
Also, |Wo,z — Go,z\ = C 5 (ir/l — x)~ 2 , for x G '[vr/2/, 6—1]. 

As before we denote 0o and <pi as Hopf differentials corresponding to 
harmonic maps Wo(t) and Wi(t), respectively. We also denote /io and fi\ as 
the corresponding Beltrami differentials in M nS and Mi PI Si, respectively. 
Write (p Go = pG 0iZ (t)G 0tZ (t) and <f) Gl = pG ltZ {t)G 1>z {t). Notice that in M fl 
Eq, map Go is the cylinder map hence harmonic, so 0g o is the Hopf differential 
corresponding to Go- When t = we have Wq = Go = identity and p — a, 
hence we can differentiate 4>q = pGo tZ (t)Go, z (t) in t at t — 0, and find that 
|0g o - <fiw 1 = o-\W . z - G 0>z \ < C 5 l 2 x~ 2 csc 2 (lx) = 0(1) for x G [a + l,n/2l], 
and \<p Go - <p Wo \ < C 5 l 2 (n/l - x)- 2 csc 2 {lx) = 0(1) for x G [n/2l,b - lj. 
Therefore we have proved 

Lemma 11. |0 Go — (j) Wo \ = 0(1) /or x G [a + 1, 6 — 1]. 

Similarly we have \<j> - <j>G \\Mir\Zi = 0(1), and \ fa - 0GillA/ o ns o = 0(1), 
also |0i - ^gJImihSi = 0(1). 
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3.5 Proof of the Main theorem 



In last subsection, we constructed families of maps and found that they are 
reasonably close to the corresponding families of harmonic maps between 
surfaces. Now we are ready to adapt the estimates in the model case to the 
general setting, and prove the proposition in this situation, which will imply 
our main theorem. 

Recall from §3.3 that the sectional curvature is represented by R/U. With 
Lemma 6 in the model case, the inequality II > iTM nE ~ still holds. 
Hence |-R|/n = 0(l 3 )\R\, so it will be sufficient to show that \R\ = o(Z~ 3 ). 
From Lemma 7, we have \R\ < 4 J E D(\ii \ 2 )\fii\ 2 adxdy, so we will show that 
Lemma 8 still holds, which implies desired curvature estimate immediately. 

Now we are about to estimate J" E D(\ii Q \ 2 )\iii\ 2 adxdy, which breaks into 
3 integrals as follows: 

/ D(\ii \ 2 )\/j,i\ 2 adxdy = / D(\^i \ 2 )\^i\ 2 adxdy 

Jt, «/M nE 

+ / D{\^i Q \ 2 )\^i\ 2 adxdy 
./MinEi 



+ / D(\fi \ 2 )\fli\ 2 adxdy (12) 

JK 



where K is the compact set disjoint from (M n S ) U (Mi Pi Si). 

For the third integral, from the previous discussion, because of the con- 
vergence of the harmonic maps to the harmonic maps of "noded" problem, 
we have both |/io| and \fii\ are bounded. The maximum principle implies 
that D(\fi \ 2 ) = H < sup{\fio\ 2 }. Note that K is compact, hence we have 
the third integral is the order of 0(1). 

From Lemma 4, we have (A — 2)H = — 2^L, so we can rewirte (12) as 

D(\[io\ 2 )\fii\ 2 o'dxdy = / H\fii\ 2 adxdy 
E iM ns ( , 



+ / n\/di\ 2 adxdy + 0(l) (13) 
J Mi nEi 



Now we will look at the first and the second integrals in (13). Recall that 

(A — 2)?i G = — 2-^2^, where 7i G is the holomorphic energy corresponding 
to the model case when the harmonic map is the cylinder map, and <f> G is 
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the quadratic differential corresponding to the constructed map G . We also 
denote jj, G to be the Beltrami differential corresponding to <f) G . We assign 
similar meanings for 0f and fif. Also recall that |0o — 00*1 — 0(1) and 
|0! - 0f | = 0(f) in (M n S ) U (Mi n Si), here 0(1) is bounded in / 
for small /. So we can set some A = 0(1) (bounded in / for small I) such 
that |0 O | 2 < A 2 |0(f | 2 and at the boundary of (M H S ) U (Mi n S x ) satisfies 

H < XH G '. For example, we can take A = l + maxd K (^, this A = 0(1) 

because at dK = 9((M nS )U(M 1 nSi)), both H, H G , and S are bounded. 
Therefore, 

(A - 2){H - XH G ) = 2 A2| ^-^° |2 > 

So (H — XH G ) is a subsolution to the equation (A — 2)Y = whose 
solutions have the form of Y(l;x) = B 3 cot(lx) + -84(1 — lxcot(lx)), where 
constants B 3 and B 4 satisfy that B 3 = 0(1) and B 4 = 0(1). Hence in 
(M n S ) U (Mi n Si), we have H < XH M + B 3 cot(lx) + B 4 (l - lxcot(lx)). 
We apply this into (13) and find that, 



/ D(\[lo\ 2 )\fii\ 2 crdxdy = / ii\i±\\ odxdy + 0(1) 

Jt, J(Af nEo)u(AfinEi) 

< / (XH G + Y(i,x))(\fi 1 \ 2 ) ( rdxdy 
i(M nSo) 

+ / (\H G + Y(l;x))(\fi 1 \ 2 )adxdy + 0(1) 

J (MiPlEi) 

< / (\H G + Y(l;x))(2\fi 1 \ 2 + 2]^- fi G \ 2 )dA 

+ / (A7i: G + F(/;a;))(2|/ i i| 2 + 2|/ i i-^| 2 )rfA 
•/ Mi nEi 



'MinEi 
+ 0(1) 



Recalling the computation in §3.3, and |/ti — /if | = |0i — <Pi\/cr, where 
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>i — <pf I < CrJ 2 x 2 csc 2 (lx) for x G [a + 1, 7r/2£], we have the following 



Af ns 



AfinSi 





/ (Aft G )|/if \ 2 adxdy = 

J Ma 


o{r 2 ) 




1 (\H G )\fif\ 2 adxdy = 

J M ± 


o(r 2 ) 




/ Y(l, x) \fi G \ 2 adxdy = 

J Mo 


0(1) 




/ Y(l, x)\fii\ 2 adxdy = 

J Mr 


o(r 2 ) 


(ft G 


+ Y(l;x))(\fn-fi G \ 2 )dA = 


o(r 2 ) 


(ft G 


+ Y(l;x))(\fn~fi G \ 2 )dA = 


o(r 2 ) 



Therefore J s D(\fiQ\ 2 )\fii\ 2 adxdy = 0(1 2 ), with Lemma 7, we have \R\ < 
4 J s D(|/io| 2 )|/ii| 2 cx(ix(i?/ = 0(l~ 2 ). Apply this into the curvature formula, we 
have |-R|/n = 0(1) — > as / — > 0, which completes the proof of proposition. 

4 Punctured Surface Case 

In previous sections, the argument on a compact surface is enough to prove 
our main theorem. However, we remark here that we allow the surface has 
finitely many punctures. 

For the case of punctured surfaces, the existence of a harmonic diffeo- 
morphism between punctured surfaces has been investigated by Wolf (|26j) 
and Lohkamp ( |14j ) . In particular, Lohkamp (|14j) showed that a home- 
omorphism between punctured surfaces is homotopic to a unique harmonic 
diffeomorphism with finite energy, and the holomorphic quadratic differential 
corresponding to the harmonic map in the homotopy class of the identity is 
a bijection between Teichmuller space of punctured surfaces and the space 
of holomorphic quadratic differentials. 

In this case, the set S\((M H S ) U (Mi n Si)) is no longer compact. 
Let Kq be a compact surface with finitely many punctures, and {K m } be 
a compact exhaustion of K . We now estimate j K -D(|/i | 2 )|/ii| 2 cL4. Let 
H(t) be the holomorphic energy corresponding to the harmonic map w(t) : 
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(K ,a) — > (K ,p(t)), then -fiT(t) is bounded above and below, and has nodal 
limit 1 near the punctures ([2H]), hence both |/i | 2 and |/ii| 2 have the order 
o(l) near the punctures. To see this, we consider K as the union of K m and 
disjoint union of finitely many punctured disks, each equipped hyperbolic 
metric z l^ 2 ^ ■ Then = 0(\z\log 2 z) — > as z tends to the puncture, since 
the quadratic differential has a pole of at most the first order. A similar result 
holds for We notice that 8Kq is the boundary of the cylinders, where 
the harmonic maps converge to a solution to the "noded" problem as I — > 0, 
hence D(\fi \ 2 ) = H(t) is bounded on dK . Therefore we apply Omori-Yau 
maximum principle ([T^] [HI]) to (A — 2)H = — 2|/i | 2 on K and obtain 
that sup(D(\fi \ 2 )) < max(sup(\fi \ 2 ) , max(D(\fi \ 2 ))\Q Ko ) = 0(1). Hence we 
have 

/ D(|/i | 2 )|/ia| 2 cL4 < / supQfiol^lfi^dA 

JK Jk 

< O(l)O(l)Vol(K ) 
= 0(1) 

In other words, our proof carries over to the punctured case, which completes 
the proof of our Main theorem. 
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